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Abstract 

We consider the stochastic heat equation with multiplicative noise u t = 
|A« + uW in R+ x R , whose solution is interpreted in the mild sense. 
The noise W is fractional in time (with Hurst index H > 1/2), and colored 
in space (with spatial covariance kernel /). When H > 1/2, the equation 
generalizes the Ito-sense equation for H — 1/2. We prove that if / is 
the Riesz kernel of order a, or the Bessel kernel of order a < d, then 
the sufficient condition for the existence of the solution is d < 2 + a (if 
H > 1/2), respectively d < 2 + a (if H = 1/2), whereas if / is the heat 
kernel or the Poisson kernel, then the equation has a solution for any 
d. We give a representation of the fc-th order moment of the solution, in 
terms of an exponential moment of the "convoluted weighted" intersection 
local time of k independent a!- dimensional Brownian motions. 
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1 Introduction 



The study of stochastic partial differential equations (s.p.d.e's) driven by a Gaus- 
sian noise which is white in time and has a non-trivial correlation structure in 
space (called "color" ) , constitutes now a classical line of research. These equa- 
tions represent an alternative to the standard s.p.d.e.'s driven by a space-time 
white noise. A first step in this direction has been made in [lOj . where the 
authors identify the necessary and sufficient condition for the existence of the 
solution of the stochastic wave equation (in spatial dimension d = 2), in the 
space of real-valued stochastic processes. The fundamental reference in this 
area is Dalang's seminal article [9], in which the author gives the necessary and 
sufficient conditions under which various s.p.d.e 's with a white-colored noise 
(e.g. the wave equation, the damped heat equation, the heat equation) have 
a process solution, in arbitrary spatial dimension. The methods used in this 
article exploit the temporal martingale structure of the noise, and cannot be 
applied when the noise is "colored" in time. Other related references are: [55] . 
[23], [33], [33] and [HJ. 

Recently, there has been a growing interest in studying s.p.d.e.'s driven by 
a Gaussian noise which has the covariance structure of the fractional Browniam 
motion (fBm) in time, combined with a white (or colored) spatial covariance 
structure. (Recall that an fBm is a centered Gaussian process (Bt)t>o with 
covariance E{B t B s ) = R H (t,s) := {t 2H + s 2H - \t - s\ 2H )/2, with H € (0,1). 
The Brownian motion is an fBm of index H — 1/2. We refer the reader to the 
expository article [27j . for a comprehensive account on the fBm.) This interest 
comes from the large number of applications of the fBm in practice. To list only 
a few examples of the appearance of fractional noises in practical situations, 
we mention [3U] for biophysics, [3J for financial time series, [T3J for electrical 
engineering, and [7] for physics. 

In the present article, we consider the stochastic heat equation with a multi- 
plicative Gaussian noise, which is fractional (or white) in time with Hurst index 
H > 1/2 (respectively H = 1/2), and has a non-trivial spatial covariance struc- 
ture given by a kernel /. As in [9J, we assume that / is the Fourier transform of 
a tempered measure \i. (Note that the particular case of a spatially white noise 
arises when / = 5q.) More precisely, we consider the following Cauchy problem: 

^ = ^Au + uW, t>0,xeR d (1) 
Uo,x = Uq{x), x e M d , 

where uo G Cf,(R d ) is non-random, and W is a formal writing for the noise 
W — {W(h);h € HV} (to be introduced rigourously in Section^]). 

Before discussing the multiplicative case, we recall briefly the known results 
related to the existence of the solution of the stochastic heat equation with 
additive noise: 

= -Au + W, t>0,xeR d (2) 
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u 0tX = o, ier, 

When H = 1/2 and f = 5q, equation © admits a solution in the space of 
real- valued processes, if and only if d = 1. This phenomenon can be explained 
intuitively by saying that, while the Laplacian smooths, the white noise roughens 
(see also [16]). If the spatial dimension d is larger than 2, then the roughness 
effect of the white noise overcomes the smoothness influence of the Laplacian. 

What happens when the space-time white noise is replaced by a noise which 
is fractional in time, but continues to be white in space? This situation has been 
studied in several papers such as [13] , [23] j [2H] > [33] and recently in 1 . In this 
case, a necessary and sufficient condition for the existence of the solution of j2]) 
is d < AH, which allows us to consider the cases d = 1, 2 or 3, for suitable values 
of H. This can be interpreting by saying that for H > 1/2, the noise roughens 
a little bit less, and the smoothness influence of the Laplacian overcomes the 
roughness of the noise. If the noise is colored in space, the conditions for the 
existence of the solution of ([2]) depend on the noise regularity in space. For 
example, if / the Riesz kernel of order a, or the Bessel kernel of order a, the 
necessary and sufficient condition for the existence of the solution of ([2]) is 
d < AH + a, whereas if / is the heat kernel or the Poisson kernel, the solution 
exists for any d > 1 and H > 1 /2 (see [T] , as well as Appendix B for a correction 
of the result of [1]). 

Another explanation of this phenomenon is given in [16j , and it is related to 
the local time of the stochastic processes associated with the differential operator 
of the s.p.d.e. In the particular case of the stochastic heat equation driven by a 
space-time white noise, the solution exists only in dimension d = 1 because this 
is the only case when the d-dimensional Brownian motion has a local time. 

We now return to the discussion of equation {!]). This equation has been 
studied recently in [18] . when the noise is fractional in time, and white in space. 
In this article, it is proved that a sufficient condition for the existence of the 
solution (in the space of square- integrable processes) is d < 2: if d = 1, then 
equation |1| has a solution in any time interval [0, T], but if d = 2, this equation 
has a solution only up to a critical point To (i.e. it has a solution in any interval 
[0,T], with T < To). It is not known if this condition is necessary as well. 
There still is a connection with the local time, in the sense that the second-order 
moment of the solution is equal to the exponential moment of the "weighted" 
intersection local time L t of two independent d-dimensional Brownian motion 
B 1 and B 2 , written formally as: 

L t := H(2H- 1) f f \r - s\ 2H - 2 S (B^. - B 2 s )drds. 
Jo J o 

In the present article, we consider equation |1| driven by the Gaussian noise 
introduced in [TJ. This noise is fractional in time with Hurst index H > 1/2, and 
colored in space, with covariance kernel / chosen among the following: the Riesz 
kernel, the Bessel kernel, the heat kernel, or the Poisson kernel (see Examples 
|2.1II2.4|) . The case of the fractional kernel f(x) = rji=i H i( 2H i ~ l)\xi\ 2Hl ~ 2 
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with 1/2 < Hi < 1 has been examined in [TJj, using methods that rely on the 
product form of /. These methods cannot be used in the present article, since 
in our case (except the heat kernel), / is not of product type. For the fractional 
kernel, it was proved in [17] that the sufficient condition for the existence of the 
solution is d < 2/(2H - 1) + £* =1 H t . 

As in the case of equation ([2]) with additive noise, we find that the existence 
of the solution depends on the roughness of the noise. If H > 1/2, and / is the 
Riesz kernel of order a, or the Bessel kernel of order a < d (which are "rough" 
kernels), then a sufficient condition for the existence of the solution is d < 2 + a: 
if d < 2 + a the solution exists in any time interval [0, T], whereas if d = 2 + a, 
the solution exists only up to a critical point To. If / is the heat or the Poisson 
kernel (which are "smooth" kernels), the solution exists in any time interval, for 
any d > 1 and H > 1/2. If / is one of the "rough" kernels mentioned above, we 
prove that if the solution exists, then d < AH + a. This shows that for H = 1/2, 
the necessary and sufficient condition for the existence of solution is d < 2 + a. 
It remains an open problem to identify the necessary and sufficient condition 
for the existence of the solution, in the case of H > 1/2. 

The existence of the solution is connected to the "convoluted weighted" 
intersection local time Lt, written formally as: 

L t = H(2H - 1) f [ f \r- s\ 2H - 2 5 Q (Bl - B 2 S - y)f(y)dydrds. 
Jo Jo Jm d 

More precisely, the second-order moment of the solution can be expressed as: 

E[u 2 J = E [u Q (x + Bl)u (x + B 2 ) exp(L t )] . 

As in [18], this expression can be extended to the moments of order k > 2, using 
k independent d-dimensional Brownian motions. 

This article is organized as follows. Section 2 contains some preliminaries 
related to analysis on Wiener spaces. In Sections 3, we discuss the existence of 
the solution. In Section 4, we examine the relationship with the "convoluted 
weighted" intersection local time. 



2 Preliminaries 

We begin by describing the kernel which gives the spatial covariance of the noise. 
As in [S], let / be the Fourier transform of a tempered distribution [i on M. d , i.e. 

fix) = [ e-**/i(dO, V.t e M d , 



where £ • x denotes the scalar product in M. d . Let V(M. d ) be the completion of 
{1a; A £ B b (R d )}, where B b (R d ) denotes the class of bounded Borel sets in R d , 
with respect to the inner product 

(1a, 1b)-p(r<<) = / / f(x-y)dydx. 
J a Jb 
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We consider some examples of kernel functions /. In what follows, |x| denotes 
the Euclidian norm of x £ M. d . 

Example 2.1 If (i(d£) — \£,\~ a dt; for some < a < d, then / is the Riesz kernel 
of order a: 

f(x)= la , d \x\- d+a , 
where la . d = T((d - a)/2)2- a ir^ 2 /T(a/2). 

Example 2.2 If = (1 + |£| 2 )~ Q/2 d£ for some a > 0, then / is the Bessel 

kernel of order a: 

poo 

f{x)=i a / w^/^e^e-^/^dw, 



where i a = (47r) Q / 2 r(a/2). In this case, T(M. d ) coincides with H~ a/2 {R d ), the 
fractional Sobolev space of order —a/2; see e.g. p. 191, [15] . 

Example 2.3 If = e -71 ""'^ 2 / 2 ^ for some a > 0, then / is the heat kernel 

of order a: 

f(x) = (2nay d / 2 e-W 2 ^ 2a \ 

Example 2.4 If n(d£) — e^^ 2 a ^ d^ for some a > 0, then / is the Poisson 
kernel of order a: 

f(x) = C d a(\x\ 2 + aT {d+1)/2 , 
where C d = 7r-( d+1 )/ 2 r((d + l)/2). 

As in [I], if H > 1/2, we let WP be the Hilbert space defined as the comple- 
tion of {l[o,i]xAi t > 0,A £ Bb(R d )} with respect to the inner product 

(l[o,t]xA, 1[o, s ]xb)hp = a H / / / / \u- v\ 2H ~ 2 f(x - y)dydxdvdu, (3) 

JO JO JAJB 

where a# = H(2H — 1). The space ftp is isomorphic to H ® -p(M d ), where ft 
is the completion of {l[o,t]j £ > 0} with respect to the inner product 

t ,s 

2H-2 r 



( 1 [o,(]> 1 [o, s ]>H="ff / / \u-v\ dvdu. 



o Jo 



If H = 1/2, we let HV be the completion of {l [0 ,t]xAit > 0,A e B b (R d )} 
with respect to the inner product 

( 1 io,t]xAi 1 io,»]xfl)«p = ( tAs ) / / f{x-y)dydx. 

JAJB 

In this case, the space WP is isomorphic to L 2 (R+) ® 7 , (K d ). 

We note that in both cases, the space TCP may contain distributions. 
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Let W = {W(h); h £ TCP} be a zero-mean Gaussian process, defined on a 
probability space (Cl,^F,P), with covariance 

E(W(h)W(g)) = (h,g)nv- 

The process W introduce formally the noise perturbing the stochastic heat 
equation. This noise is considered to be "colored" in space, with the color given 
by the kernel /. If H > 1/2, the noise is fractional in time, whereas if H = 1/2 
the noise is white in time. 

We now introduce the basic elements of analysis on Wiener spaces, which 
are needed in the sequel. For a comprehensive account on this subject, we refer 
the reader to [SB] and [2"5] , 

We begin with a brief description of the multiple Wiener (or Wiener-Ito) 
integral with respect to W. Let T w be the cr-field generated by {W(h);h £ 
TCP}, H n (x) be the n-th order Hermite polynomial, and TCP n be the closed 
linear span of {H n (W(h))\ h £ TCP} in L 2 (tt, T w , P). The space TCP n is called 
the n-th Wiener chaos of W. 

It is known that L 2 (VL, F W ,P) = ®%L HV n , and hence every F £ L 2 (VL, T w , P) 
admits the following Wiener chaos expansion: 

oo 

F = Y / MF), (4) 

where J n : L 2 (n, F w , P) — > TCP n is the orthogonal projection. By convention, 
TLVo = K and J (F) = E(F). 

For each n > 1, and for each h £ TCP with ||/i||-k-p = 1, we define 

I n (h® n )=n\ H n (W(h)). 

By polarization, we extend /„ to elements of the form hi ® . . . ® h n (see p. 230 of 
[T7] : e.g. h 1 ®h 2 = [(hi + h 2 )® 2 - (hi - h 2 )® 2 ]/^)- By linearity and continuity, 
we extend the definition of I n to the space TCP® n . (Note that if \e-i\i > 1} is 
a CONS in HP, then {e^ ® . . . ® e in ; ij > 1} is a CONS in W® n .) For any 
h £ TCP® n , we say that 

I n (h) ■= / h(ti,Xi,...,t n ,x n )dW tl , xl ...dWt n} x n 

J(R + xK d )™ 

is the multiple Wiener integral of h with respect to W. We have 

E(I n (h)I n (g)) = n! (h,g) nv ^, Vh,g£ TCP m , 

where fo(*i,a;i, . . .,*„,£„) = (n!) _1 Eaes„ ^(^(l) , ^a(l) , • ■ ■ ,t<r(n), ^(n)) is the 
symmctrization of ft. with respect to the n variables (ti, Xi), . . . , (t n , x n ), and S n 
is the set of all permutations of {1, ... ,n}. By convention, we set Io(x) = x. 

The map I n : TCP® n — > is surjective. Moreover, for any F n £ TLP n , 

there exists a unique /„ € TCP® n symmetric, such that I n (f n ) — F n . Using (TJ|), 
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we conclude that any F £ L 2 (Q, T % P) can be written as: 

OO 

F = £/„(/„), (5) 



where /o = E{F) and /„ G HV® n [ s symmetric and uniquely determined by F. 
We have: 

OO OO 

E\F\* = J2m(fn)\ 2 = J2 nl WfnWnvvn- 

n=0 n=0 

We now introduce the stochastic integral with respect to W. Let u = 
{ut, x \ (t,x) G M+ x M d } be an JF^-measurable square-integrable process. By 
©, for any (t,x) el + x M d , we have 



E(u t , x ) + J2l n (f n (;t,x)), (6) 



where f n (-,t,x) G TCP® n is symmetric and uniquely determined by Ut jX - 

For 

each n > 1, let /„ be the symmetrization of f n with respect to all n + 1 variables. 
Let /o = E(u). We say that u is integrable with respect to W if /„ G 
np ®(n+i) for every „ > o, and E^=o J n+i(/n) converges in L 2 (f7). In this 
case, we define the stochastic integral 



/•OO v 

8(u) := / u s <5W 5 = V/„+i(/n) 
J o „ =0 



Note that: 

OO 

E\S(u)\ 2 = J2(n + iy- \\fn\\ 2 Hv ^ + v 

The following alternative characterization of the operator 6 is needed in 
the present article. Let S = {F = f(W(hi),...,W(h n ));f G C^(M. n ),hi G 
7iV, n > 1} be the space of all "smooth cylindrical" random variables, where 
C£°(R d ) denotes the class of all bounded infinitely differentiablc functions on 
M. n , whose partial derivatives are also bounded. The Malliavin derivative of an 
element F — f(W(hi), . . . , W{h n )) G S, with respect to W, is defined by: 

DF := J2 J^^ 1 )' ■ ' • ' W ( h n))hi- 

l — l 1 

Note that DF G L 2 (tt; HP); by abuse of notation, we write DF = {D t . x F; (t, x) G 
[0,T] x M d } even if D t , x F is not a function in (t, x). We endow S with the 
norm |jF|j 2 1-2 := E\F\ 2 + J5||DF||^ P) we let ED 1 ' 2 be the completion of 5 
with respect to this norm. The operator D can be extended to ED 1 ' 2 . Then 
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S : Dom S C L 2 (fl;H'P) — > L 2 (f2) is the adjoint of the operator £>, and is 
uniquely defined by the following duality relationship: u G Dom 6 if and only if 

E(FS(u)) = E{DF, u) HV , VF e B 1 ' 2 . (7) 

Note that u e Dom 6 if and only if u is integrable with respect to W . (In the 
literature, 5 is called the Skorohod integral with respect to W.) 

3 Existence of the solution 

In this section, we give conditions for the existence of the solution of equation 
Let pt(x) be the heat kernel on R d , i.e. 

^) = ^-p(-^), t >0,*eR- 

For any bounded Borel function ip : W d — > M, let pt<f{x) = j^ d pt(x — y)ip(y)dy. 
For each f > 0, let JF t be the er-field generated by {W(l[o a ]x.A)j s £ [0,t],ie 

The solution of equation (TTJ) is interpreted in the mild (or evolution) sense, 
using the stochastic integral introduced above. More precisely, we have the 
following definition. 

Definition 3.1 An {Tt)t- adapted square-integrable process u — {tit iX ; (t,x) £ 
M + x R d } is a solution to (QP if for any (t,x) S M+ x R d , the process 

{Yl'y = l m (s)p t - s (x - y)u SlV ; (s,y) e M+ x M d } 

is integrable with respect to W , and 

u t , x = PtM x ) + [ [ Y sySW SyV . 
Jo Js. d 

By (|7|), the above definition is equivalent to saying that for any (t, x) € K+ x M. d , 
Ut,x G L 2 (ty, u t,x is ^(-measurable, and 

E{u t<x F) = E{F) Pt u (x) + E{Y*> X , DF) HV , VF e B 1,2 . (8) 

The next result establishes the existence of the solution u = {u t , x \ (t,x) € 
R + x R d }, as a collection of random variables in L 2 (£l). As in [T5] (see also [5J, 
EH) [25] , [30], [32], [35]), one can find a closed formula for the kernels f n (-,t,x) 
appearing in the Wiener chaos expansion ([6]) of u t , x - 

Proposition 3.2 In order that equation (Q]) possesses a solution it is necessary 
and sufficient that for any (t, x) € R+ x R d , we have 

oo 

J^n! || /«(•,*, a) |Ihp®*> < oo, (9) 

n=0 
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where 

1 - 

/„(*!, Xi, . X„,t,x) = — l[Pt plj+1) -t Kj) (Xp(j+l) - X pU) )p tp(1) Uo(x p(1) ), 



p denotes the permutation of {1, 2, . . . , n} such that t p m < tp(2) < ■ ■ ■ < t p ( n ), 
tp(n+i) = t and Xp( n +\) = x. In this case, the solution u is unique in L 2 (fl), 
admits the Wiener chaos decomposition and 

oo 

£K*I 2 = X>I \\U;t,x)f HV ® n . ( 10 ) 

ra=0 

We begin to examine condition Note that 

a n (t,x) := (n!) 2 ||/„(-,t,a:)||^ 9 - < hof^a^t) , (11) 
with equality if Uq = 1, and hence 



oo ^ oo ^ 

E\u t , x \ 2 = V —.a n {t, x) < ||mo||L V —.a n (t), 

where 



„ n! * — ' n\ 

n=0 n=0 



^/i„, (F »n; =1 \ S] -t 3 \ 2H - 2 r {n) (s,t)d S dt if h > 1/2 

a n {t) = { (12) 
/ [0it]B r (n) (s,s)ds if if =1/2 

and 

V>* (n) (s,t: := / [[ <!i)llPt pU+1) -t, u) (x rlr - lr v> - r ril ,) 

3= ' 



s,t) := / n^~ y ^n^ 

• /R2 " d .7 = 1 .7 = 1 



II^u+d-'Mj) (^(i+i) - y*(j))d*dy. (13) 
i=i 

In the above integrals, we denoted s = (s%, . . . , s n ), t = (ti, . . . , t n ),x — (x%, . . . , x n ) 
y = (j/i, . . . , y n ), and we chose the permutations p and a of {1, . . . , n} such that 

< tp(i) < t p (2) < ■ ■ ■ < tp( n ) and < s .( 1 ) < s ff ( 2 ) < . . . < s a{n) , (14) 

with <p(„+i) = s CT („+i) = i and x p („+i) = y a (n+i) = x. 
Note that 

V>* (n) M) = (gi n \gi n) ) vm ®n, Vt.se [o,t]» 

where 

n 

^ n) (a;i,...!C n ) = Y[p tK . +1) _ tp{j) (x p{j+1) -x p{j) ) 
4=1 
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and the permutations p and a are chosen such that (fl"4")) holds. 
As in pQ, for any y, 2 G M. d and u, w > 0, we denote 



Jf{u,v,y, z) := / / p u {x - y)p v (x' - z)f(x ~ x')dxdx' . 

Lemma 3.3 (i) If f is the Riesz kernel of order a, or the Bessel kernel of order 
a < d, then 

J f (u, v, y, z) < D a , d (u + vY {d - a ^ 2 , Vy, z E R d , 

where D a d is a positive constant depending on a and d. 

(ii) If f is the heat kernel of order a, or the Poisson kernel of order a, then 

Jf(u,v,y, z) < C a ,d, Vy, z e R d . 

Proof: Note that Jf(u,v,y, z) — E[f(y — z + y/uY — i/vZ)], where Y and Z 
are independent <i-dimensional standard normal random vectors. We use the 
following inequality: (see (3.24) of pQ) 

B [ e -|»-,+>/Ky-V3Uz|V(4»)] <( 1 + ^ . (15) 



(i) In the case of the Riesz kernel, this inequality has been shown in the 
proof of Theorem 3.13 of pQ. Suppose now that / is the Bessel kernel of order 
a < d. Using iflgl). 



J f (u,v,y,z) = i a / w ^- d )/^e~ w E[e^- z+ ^ Y -^ z \ 2 '^]dw 
Jo 

< y a J™ w (°- d )/^ e - w (^ + ^-^ 2 dw 

- < f wa/2 - le ' w (» + u -¥) ~ m dw - ^ (T 

where I a ^(x) := J °° w a / 2 ~ l e~ w (w + x)~ d l 2 dw. The result follows, since 
I a ,d( x ) < x- d/2 w a/2 ~ 1 e- w dw+ w a/2 - 1 e- w (w + x)- d/2 dw 

Jo Jx 

f 1 /"OO 

= x- d ' 2 x a ' 2 J y a ' 2 - 1 e- x ydy + x-^' 2 J y al2 ~ x e~ xy {y + iy d/ 
< x-^/ 2 / y^^dy + x-^-^ 2 / y-^/^dy = K a . d x^ d 



JO Jl 

where K a ^ — 2/a + 2/(d—a) = 2d/[a(d— a)], and we used the fact that a < d. 
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(ii) If / is the heat kernel, using P^|) . we obtain: 



J f (u,v,y,z) = (2™)- d / 2 £[ e Hy-2+^-^l 2 /(2a)] 

< (2TO )-Wi + !i±£)~' ,/a 

= (27r)- d / 2 (a + u + v)- d/2 < (2ira)- d / 2 . 

If / is the Poisson kernel, we have: 

J f (u,v,y,z) = C d aE[(\y-z + V^Y-V^Z\ 2 + a 2 )-^ +1 ^ 2 } 
< C d a(a 2 r id+1)/2 = C d a- d . 

□ 

Lemma 3.4 (i) If f is the Riesz kernel of order a, or the Bessel kernel of order 
a < d, then for any s, t S [0, t] n , 

V>* (n) (s,t) < (D a4 2-^ d - a ^ 2 Y [l3(s)[3{t)Y (d ~ a)/ \ 

where (3{s) := Y[ n j=1 {s a(j+1) ~ s a(j)), P(t) '•= n"=i(*p(i+i) ~ *p(j))> and the 
permutations p and a are chosen such that \1J$ holds. 

(ii) If f is the heat kernel of order a, or the Poisson kernel of order a, then 
for any s, t £ [0, t] n , 

^ (n) M)<C" d , 
where C a ^ d is a constant depending on a and d. 

Proof: By the Cauchy-Schwartz inequality, 

^* (n) (s,t) < ^*< n >(s,s) 1 /V (n) (M) 1/a - 

To find an upper bound for t/»*w (s, s), we use Lemma 13.31 to estimate the 
following integrals: 



Pu 3 {xa(j+i) ~ x a{j ))p Uj (y a{j+1) - y a (j))f{xa(j) - Va^dx^dy^ 

i 

= Jf( u j> u j> x a(j+l),ya(j+l))> w ith Uj = s <r(j+l) - s a(j), j = 1> • • • , «• 

(i) In this case, l i < D a . d [2{s <j{:j+l) - s CT (j))]~ (d ~ a)/2 and 

-{d-a)/2 



r {n \s,s)<Dl d 2 



-n(d-a)/2 



3 = 1 



(ii) In this case, I 3 < C aA and ip* 'W (s, s) < C^ d - □ 
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If H > 1/2, it was proved in [55] that there exists f3jj > such that 

/>00 /• X / / x \ '2-1 f 

10 JO 

for any <p £ L 1 / H (R + ). Hence, 



tp(8)tp(t)\t - s\ 2H - 2 dsdt <$bU \<P(t)\ 1/H dt) 



a n H J^ ^(s)^(t) f[ \tj - Sj \ 2H - 2 dsdt < 0* (J Mt^dt j 



2H 



(16) 



3=1 

for any 93 € i 1 /^ (M™). If i? = 1/2, we let (3 H = 1. 
We need the following auxiliary result. 

Lemma 3.5 Lei 

/„(*, ft) := / [(* - s n )(s n - s n _i) . . . (s 2 - si)]' l ds, 

J T n 

where T n — {s = (si, . . . , s n ); < si < s 2 < ... < s„ < <}. Tften /„(i, ft) < 00 
if and only if 1 + ft > 0. In iftis case, 

T (th)= r ( 1 + fe ) n+1 f n(l+/t) 

r(n(l + ft) + l)* 

Proof: First note that J S2 (s 2 - si) h dsi = s% +1 /(h + 1), and then 

/ 3 (s 3 ~ s 2 )' l S2 +1 ds 2 = sf + 2 /3 (ft + 2, ft + 1) 
io 

= sl {h+1) f3((h+l) + l,h + l)) 

where P(a,b) = J Q a; a_1 (l — x) b ~ 1 dx is the beta function, and we used the 
change of variables s 2 /s3 = z. In this way, the integral I n {t, ft) becomes 

I n (t,h) = j—/3((h + l) + l,h + l)...(3((n-2)(h + l) + l,h + l) 
' ' s (n-iKh + i) {t _ Sn) h dSn 

= t r ^ h+ ^l3 ((ft + 1) + 1, ft + 1) . . .0 ((n - l)(ft + 1) + 1, ft + 1) . 

Using the fact that (3(a,b) = T(a)T(b)/T(a + b) for a,b > and T(z+1) = zT(z) 
for any z > 0, we obtain the desired conclusion. □ 

Using Lemma EP1 Lemma [331 and we obtain the following estimate for 

Qn(t)- 

Proposition 3.6 Suppose that H > 1/2 and Zet a„(t) 6e jroen 6y i f 
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(i) If f is the Riesz kernel of order a, or the Bessel kernel of order a < d, 

and 

H>^, (17) 

then 

a n (t) < C* H ^ a D{t) n {n\f d -^/\ for anyt>0,n> 1, 
where C* H d > is a constant depending on H, d, a, and 

/ , x 2H , j x -[2H-(d-a)/2] 

D(t) = D a , d 2-(«-<*» 2 f H T (l - j (l - fH- (d - a) /2_ 

(ii) If f is the heat kernel of order a, or the Poisson kernel of order a, then 

a n (t) < C(t) n , for any t > and n > 1, 
where C{t) = C a4 t 2H . 

Remark 3.7 Proposition 13.61 shows that X)n a «W/ n - grows exponentially in 
time in some cases, and faster than exponentially in other cases. 



Proof of Proposition [3761 We only give the proof in the case H > 1/2, the 
case H = 1/2 being similar. We use the definition ([12"]) of a n (t). 

(i) Let h = -(d- a)/(4H). By Lemma I3H(T) and dTBD , we obtain: 



P n 

,(*) < fo M 2-( d -«)/ 2 ) a" ff ]l\t 3 - s 3 \ 2H - 2 [p(s)l3(t)r (d - a)/ *dsdt 

( \ 2H 

= [D a!d 2-( d -^ 2 f3 2 H ) n (ni) 2H I n (t,h) 2H . 



Using Lemma 13.51 we obtain: 



*»(*) < r(i + hf H (D a , d 2-^)/ 2 f3 2 H ) n (nl) 2H { r(^|^ 1} t" (1+h) 
= r(i + ^) 2ff {^ Q , d 2-( d - Q )/ 2 /3|p (i + /o 2ff < 2 «d+' 1 ) } 



2ff 



IB 



r(n(l + A) + l), 



The result follows by using (3.19) of [18] , 
(ii) By Lemma [3T4J (ii) , 



a n (t) < C™ d a n H 

□ 



/ II l*i - t,\ 2H - 2 dsdt = ci d t 2Hn = c(ty 



Using Proposition 13. 6i we examine the existence of the solution of equation 
([1]). The next result is an extension of Proposition 4.3 of [18] to the case of a 
colored noise W. 
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Proposition 3.8 (i) Let f be the Riesz kernel of order a, or the Bessel kernel 
of order a < d. Suppose that either 



if > 1/2 and d<2 + a, (18) 



or 



# = 1/2 and d<2 + a. (19) 
Then fip has a unique solution in [0, T] x M. d , provided that T < Tq where 

-1/(2H-1) 

ifd<2 + a 



T =! {(l-i)V*(l-i) M [ 
[ oo 



ifd = 2 + a (2Q) 



(ii) Let if > 1/2, and f be the heat kernel of order a, or the Poisson kernel 
of order a. Then flp has a unique solution in R + x M. d . 

Remark 3.9 Either one of conditions (JTSJ) or (|T9^) is stronger that (JT7J) . 

Remark 3.10 Proposition 13.81 shows that in the case if = 1/2, the dimension 
d = 2 + a cannot be attained. 

Proof of Proposition I3T81 We apply Proposition 13. 2\ using Proposition 13.61 

(i) We have: 

^n\\\f n {-,t,x)f HV ^ < \\u \\ 2 x J2ni an ^ ~ ^™ C *H^dJ2 (n jU-(d-a)/2 - 

If d — a — 2, then the last sum is finite if D(t) < 1, which is equivalent to 
saying that t < Tq. If d — a < 2, then the last sum is finite for any t > 0, by 
Stirling's formula and D'Alembert criterion. 

(ii) We have: 

oo oo ^ oo (JU\n 

5>! \\U;t,x)\\ 2 nv0n < \\uoWl, ~^n{t) < hollLE^T <(X) - 

n=0 n=0 n—0 

□ 

The next result shows that (TTTJ) is a necessary condition for the existence of 
the solution. 

Proposition 3.11 Suppose that if > 1/2 and f is either the Riesz kernel or 
order a, or the Bessel kernel of order a. If equation (QP with uq — 1 has a 
solution in M + x R d , then (T7\ ) holds. 

Proof: Note that E\ut, x \ 2 = Y^Lq ct n (t)/n\ < oo implies that ai(t) < oo, 
which in turn implies (jTTJ) (see Appendix A). □ 
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Remark 3.12 Proposition ^. 81 and Proposition 13 . 1 II show that, if H — 1/2 and 
/ is the Riesz kernel of order a, or the Bessel kernel of order a < d, then the 
condition d < a + 2 is necessary and sufficient for the existence of the solution 
of ([T]). It remains an open problem to see if this condition is necessary, when 
H > 1/2. To resolve this issue, one needs to develop a full analysis of the range 
of a n (t) , which would include the identification of suitable lower bounds. Such 
analysis will be the subject of future investigations. 

Remark 3.13 The case H < | also constitutes an interesting line of investi- 
gation, which will be pursued in a subsequent article. We mention that in this 
case even the stochastic heat equation with linear additive fractional-colored 
noise has not been solved. The technical difficulties that appear here are re- 
lated to the structure of the space HV and the lack of the expression of the 
scalar product in this space as (|3|). Indeed, when H < h, assuming that the 
noise W(t, x) is defined for t E [0,T] and x E K d , the space TCP can be de- 
scribed as the space of measurable functions ip(s, x), s E [0, T], x E M d such that 
K*cp E L 2 ([0, T}) ® V(R d ), where 

,T ,T 

K*ip(s, x) = K{T, s)ip(s, x) + / / ((ip(r 7 x) — ip(s,x)) d\K(r, s)dr. 



So, it is necessary to use the transfer operator K* and to check (in the case of 
the additive noise) that K*g t . x E L 2 ([0, T]) <g> V(R d ) where g t ,x{s,y) =p t - s (x- 
y)l[0,r] which is in principle rather technical (in the of the stochastic heat 
equation with multiplicative fractional-colored noise, one needs to deal with the 
tensor product operator (if*)®" which has a complicated expression). 

Remark 3.14 If H = 1/2 and / is an arbitrary kernel, it was proved in |9J 
(using different methods) that the sufficient condition for the existence of the 
solution in K + x M. d of (JTJ) with vanishing initial conditions (i.e. uq = 0), is: 

^ <oo. (21) 



1 + iei 2 

(see Remark 14 of [9 ). When / is the Riesz kernel of order a, or the Bessel 
kernel of order a, (|2ip holds if and only if d < a + 2. Combining this with 
Remark 13.121 we conclude that, in the case of the two kernels, condition J21j) 
is also necessary for the existence of the solution. For an arbitrary kernel /, 
it is not known if (|2"Tj) remains a necessary condition for the existence of the 
solution. 

If H > 1/2 and / is the Riesz kernel of order a, or the Bessel kernel of order 
a, the necessary and sufficient for the existence of the solution of the stochastic 
heat equation with linear additive noise is d < AH + a, whereas if / is heat or 
the Poisson kernel, this equation has a solution for any d (see [1] and [2]). 
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4 Relationship with the Local Time 

In this section, we identify a random variable L t , denned formally as a "convo- 
luted intersection local time" of two independent d-dimensional standard Brow- 
nian motions, such that 

a n (t) = E(L«), Vn>l. (22) 

An immediate consequence of (fT0|) , (jlip and (|22|) is that the second moment 
of u tiX is bounded by the exponential moment of L t : 

£K,| 2 < ||«o||So E = W u oW 2 oo E ^y^") = ll«o||L^(e Lt ), 

71 = ' 71=0 

with equality if uq = 1 . 

To show (|22"T) . we approximate a n (t) by {a„ !6 (i)} e >o, when e — > 0, where 
the constants a n ,e(t) are chosen such that: 

a n>e (t) = E(L 7 t l e ), Vn > 1, 

for a certain random variable Lt, £ - 

To identify the approximation constants a n , e (t), we recall the definition (|12p . 
which says that a n (i) is the weighted integral of the function (s,t). The 
next lemma gives the exact calculation for the integrand ip*( n '(s,t). 



Lemma 4.1 We have 

r in) (s,t) = (2 



n)- nd I exp J -1 V ^ • & i M(d£i ) ■ • • MdCn), 



where a* k :— (t — Sj) hit — s&) + (i — hit — tk)- 



Remark 4.2 Lemma 14.11 gives a generalization -and a minor correction- to 
Lemma 4.2 of [18]. The correction refers to the fact that the result of [18] is 
stated incorrectly with the constant ajk = Sj A Sk + tj A tk, instead of a*j k . 
However, a trivial change of variables s'- := t — Sj,ij :=t — tj in the definition 
(fl"2"|) of a n (t) shows that this minor error does not affect the calculation of a n (t). 
We have indeed: 

f a&/ t o jt ] a -Il?=il*i-*il aff -V n) (s,t)dsdt if > 1/2 

«„(*) = < (23) 

[ / [M „V (n) (s,s)ds if fl" = 1/2 

with 

V> (n) (s,t) = (27r)-" d / ,x P ! - \ • «7 Jfc & • a > M(^l) • •-M(^n) 



f 1 " 1 

S ~2 E °>' fc & ' ( 
{ j,k=i J 



iP< n) {tl -s,tl - t), where 1 =(!,...,!) G M". 
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Proof of Lemma I4.lt Note that 

where T denotes the Fourier transform. Hence, 



V* (n) (M) = (</<»>, <#°>7>(r*)« 



— nd 



(27T)- 



» 



It was shown in the proof of Lemma 4.2 of [TH] that: 

-i^.[ x -(Bl-Blj] 



rg.(ti,...,£ n ) = E 



if 

3=1 
n 

if 

3=1 



-i^i-N-(B t 2 -B t 2 J] 



where B 1 = (Bj)t>o and B 2 = (B 2 ) t >o are independent d-dimensional standard 
Brownian motions. Hence, 



r {n \s,t) = (2tt)- 



IF 



«i-[(Bi.-B« 1 )-(B?.-B?)] 



ju(dfi) . ..(J,(d£ n ) 



We begin to evaluate the integrand of the above integral. We denote £j = 
• ■ ■ , &,<0, B* 1 = (B} A , . . . , Si,,) and Sf = (S^, . . . , B 2 d ). We observe that 
for any j = 1, . . . , n fixed, the random variables 

(S^ - Si ; ) - (S* f , - B*,), I = 1, . . . , d, are i.i.d, 

with the same distribution as (6,. — 6j) — (6f. — &f), where 6 1 = (b\) t >v and 
6 2 = (6 2 )t>o are independent 1-dimensional standard Brownian motions. Hence, 



IF 

3=1 



-i^-[(K- B D-( B l- B t)] 



if 



if 

3=1 



-^•,l[(B],., ! -Bt, 1 )-(B?,-, ! -B t 2 , ! )] 



if 

3=1 



^,i[(6a,-6j)-«.-6 2 )] 



if 



I I ex P { ~ 2 Y. a jk(3,ltk,l \ 

i=i y i,k=i ) 
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where for the second last equality we used the fact that the vector 

((£ - b\) - (b 2 ti - b 2 ), (bl - b\) - (bl - b 2 )) 

has a normal distribution with mean zero and covariance matrix {<J* k )\<j,k<n- 
This concludes the proof of the lemma. □ 

In what follows, we use the alternative definition (f2"5)) of a n (t), given in 
Remark 14.21 The idea is to find a suitable approximation for the integrand 
ip( n ' (s,t), by replacing the Dirac function 5q(x) with the heat kernel p e (x). 
This approximation turns out to be: 

^ n >(B,t) :=E [ p £ « - Bl - yi) . ..p E « - B 2 - y n )f( yi ) . . . f(y n )dy 

where B 1 = (B}) t >o and B 2 — (B 2 ) t >o are independent d-dimensional standard 
Brownian motions, and we denote y = (yi, . . . , y n ). 
More precisely, we have the following result. 

Lemma 4.3 Suppose that /i(<i£) = g{^)dS,, i.e. f = Tg. Then 

/[ , n n 1 

ex p \ - 2 E • & - f E i^i 2 r ■ • ■ ^ 

^ d) " [ J,A=1 3=1 J 

Proof: We first calculate the inverse Fourier transform of p £ * f: 

T-\Ps * /)(£ = F^PsiO F-'fiO = {2n)- d e-^ 2 ' 2 g{0- 
This shows that (p s * /)(o;) = (2n)- d T[e- E ^ 2 ^g{Q](x), i.e. 

|K - y|2/(2£) /(y)^ = (2n)- d [ e-^e-^ 2 / 2 g(0dC (24) 



(2tt£)<V 2 

Using ([231) with a: = B\. - B 2 . , we obtain: 



1 -IB,, -S?. -J/j| 2 /(2e) 

(2^j^ 6 ■"■ ,/ "'"" 



Therefore, 
^")(s,t) = E 



n 

3=1 



(27T£) d / 



(27r)-" d £; 
(27r)-" d 



n 

3=1 



e -^<Bl r B^) e _ el ^ /2g{ ^ 



E 



>d)n 



3=1 



-itj-(Bl-Bl) 



(2tt) 



- nc£ 



M(^l) • --K d £n)- 
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□ 



Remark 4.4 Note that the function h — p £ * / is continuous. To see this, let 
(x n ) n C M d such that x n — * x. By (f2"4"|) and the dominated convergence theorem, 
it follows that h{x n ) -> To justify this, note that \e~^' Xn e~ e \^ 2 ^ g(^)\ < 

e" e|C|2/2 5(C) for any £ <E M d ,7i > 1 and 

(2^r 2d / e- £| « |2/2 5(0^ = / |^P e /a(0l a fl(0de = Ibs^lll^) < oo, 

JR d JR d 

since p e/2 E L 2 (R d ) C 

We are now ready to define the approximation constants a n>e (t): 

( «&/ Pjt]a -n?=ii*i-*ii 2K -M B) (s,t)dBdt ifjj>i/2 

[ / [M „# ) (8,s)dB if jy = 1/2 

Note that a„ lE (i) = E(L" E ), where 

f « ff /„* r 4 J Rd |r - S | 2 ^- 2 p £ (^ - B s 2 - y)f(y)dydrds if if > 1/2 

1 Jo /h- P^Bl - B 2 - y)f(y)dyds if H = 1/2 

The random variable L tj£ is an approximation of the "convoluted intersection 
local time" L t , written formally as: 




a H Si Jo / R - k - s| 2ff - 2 <5o(^ - S 2 - y)f(y)dydrds if fl" > 1/2 
Jo /r- ^(B. 1 - Bl - y)f(y)dyds if = 1/2 



Remark 4.5 We mention that this approximation procedure has been inten- 
sively used in several papers dealing with the chaos expansion of the local time 
and Tanaka's formulas for Brownian motion (see e.g. [31 j ) or fractional Brow- 
nian motion (see [8]). Recall that the local time of the Brownian motion can 
be formally written as L(t,x) = f Sq(B s — x)ds where Sq is the delta Dirac 
function. Usually, to obtain the chaos expansion of L(t,x) one approximates 
Sq(B s — x) by the Gaussian kernel p e . 

More generally (and for the sake of a result encountered later in the sequel) , 
if T) : [0, t] 2 — > R + is an arbitrary function such that r](r, s) = rj(t — r,t — s) for 
all r, s € [0, t], we define 




rj(r, s)p £ {Bl — B 2 S — y)f(y)dydrds. 
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Thena„, e (t,7y) = E(L t , £ {rj) n ), where a n , s (*, ??) := /rotia«.rii=i r ?( s j>*j)V'e (s,t)dsdt. 



;[o,t]=" iij 
Let 

« n 

a n (t,r]):= TT r?( Sj , ^ (n) (s,t)(fedt, (25) 

■/[0,i] 2 "^l 

and note that Q! niE (t, 77) < a n (t, rj) for all e > 0. Note that a n>£ (t, rj) < a n>£ >(t, rj) 
if < e' < e. 

Lemma 4.6 Let t > be arbitrary, a) If ct2{t,rj) < 00, then 

lim E{L t<£ (rj)L tj5 (rj)) = a 2 (t, rj), (26) 

£,510 

and there exists a random variable Lt(rj) := lim £ j,o Lt, e (v) * n L 2 (fl). 

b) If a n (t,rj) < 00 for all n > 1, then the random variable L t , defined in 
part a), is p-integrable for any p > 2, and 

limE\L t , e (r 1 )~L t (r 1 )\ p = 0, forallp>2. (27) 

ej.0 

In particular, E(L t (rj) n ) = lim e ^o E(L tt£ (r)) n ) = a n (t, rj) for all n > 1. 

The random variable Lt(rj) defined in Lemma [4.61 depends on B 1 and B 2 , 
and could be denoted by: 

Lf' B \v)=[ [ [ r){r,s)8 Q {Bl - B 2 S - y)f{y)dydrds. 

JO JO JR d 

This notation emphasizes dependence on B 1 , B 2 , and the formal interpretation 
of L t (rj) as a "convoluted intersection local time" of B 1 and B 2 . 

Proof of Lemma 14. 6t As in [TB] , the proof follows by classical methods. We 
include it for the sake of completeness. To simplify the writing, we omit rj in 
the arguments below. 

a) Note that E(L i>e L t ,s) = L i4 i4 h)v( s 2, h)ipf} (s, t)dsdt, where 



4 2) s (s,t) := E 
= (2' 



f Pe « - B\ - yi )p E (Bl 2 - B 2 2 - V 2)f{vi)f{V2)d yi dy: 
V) 2 [ 2^ 2 2 j 



proof of Lemma [OJ) Then lim e , 5i0 V^fi *) = V> (2) ( s ,t)- Relation (T5SD follows 

by the dominated convergence theorem, since ipfl (s, t) < i// 2 " 1 (s, t) for all e, S > 
0, and 

rj(si,t 1 )n(s 2 ,t 2 )ip^ (s,t)dsdt = a 2 (i) < 00. 

[o,*] 4 
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^From here, we also infer that lim e |Q E(L^ E ) = ct2(t), and hence 
lim E\Lt* - L t ,s\ 2 = UmB(L? , B ) + IjmE(L? i5 ) - 2 lim .E{L t , e L t , s ) = 0. (28) 

Let (e„) n | be arbitrary. From (|28[) . it follows that (L t , £n ) n is a Cauchy 
sequence in L 2 (S1). Hence, there exists L t € L 2 (f2) such that £J|i t)En — i t | 2 — * 0. 
If (e' n ) n | is another sequence and E\L t . £ > — L' t \ 2 — ► for some L' t S L 2 (J7), 
then E\L t - L' t \ 2 < E\L t - L t . £n | 2 + £|£ Mn "- Lt.e'J 2 + E\L t< - L' t \ 2 -» 0, i.e. 
B|it — i' t | 2 = 0. This shows that L t does not depend on (s n ) n . 

b) Let p > 2 be fixed. Let (e n ) n I be arbitrary. We will prove that 
E\Lt t£n — L t \ p — > 0, by using the fact that, in a metric space, x„ — > a; if and 
only if for any subsequence N' C N there exists a sub-subsequence A^" C N' 
such that x„ — > x, as n — > oo, n € A^" (see e.g. p. 15 of [1]). 

Let AT' C N be an arbitrary subsequence. By part a), as n — > oo, n € AT', 
Li, gn — > it in L 2 (f2). Hence, Lt,e n — > it in probability, and there exists a sub- 
subsequence AT" c AT' such that Lt, En — ► it &-s., asm oo, n £ Af". Note that 
(Li, e ) e >o is uniformly integrable, since 

sup£ , (L™ e ) = supa n . E (t) < a n {t) < oo, for n>2. 

e>0 ' e>0 

By Theorem 16.14 of 5 , it follows that |L f | p is integrable and E\L t , E „ -L t \ p -> 0, 
as n — > oo , n € AT" . □ 

The next two results are the analogues of Propositions 3.1 and 3.2 of [TB] in 
the case of a colored noise. We denote $(x, a) = X^^Lo £™/( n D a f° r £ > and 
a > 0. Note that $(x, a) < oo if and only if a > 0, a; > or a = 0, x £ (0, 1). 

Proposition 4.7 Suppose thatn : [0,i] 2 — > R+ satisfies the following condition: 
\\r]\\i,t '•= max I sup / r](r,s)dr i sup / ?y(r,s)(is] < oo. (29) 

y«6[0,t]Jo re[0,t]Jo / 

// / is i/ie Riesz kernel of order a, or the Bessel kernel of order a < d, 
and d < 2 + a, then lim e m L tE {rf) = L t (rj) exists in L P (Q) for all p > 2, and 

supE [cxp (XL t ,e(v))} < C* d $ (\D(t), 1 - ^— - j , for all A > 0, 
where C* d is a constant depending on a and d, and 

/ A \ / A \ -^-(d-a)/ 2 } 

D(t) = D a , d 2-( d -^ 2 \\ V \\ u T (l (l - *!-(--)/». 

fn) /// is i/ie /ieai kernel of order a, or the Poisson kernel of order a, then 
lim e ^o Lt i£ (rj) — Lt(rj) exists in L P (Q), for all p > 2, and 

sup E [exp (\L t , e (rj))] < exp(AC(i)), /or a// A > 0, 

e>0 

wftere C(i) = C Q ,d|M|i,t^ 
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Proof: We use the definition (j2"5|) of a n (t, rf), and Lemma [3TH for the estimation 
of ^*W(s,t). 

(i) Let h = — (d — a)/2. Using the Cauchy-Schwartz inequality, condition 
(|2T)|) , and Lemma 13. 51 we get: 



r n 

v y -/[o,t] 2 " J= i 

/» n 

< (r>a,d2- (d - Q)/2 ) n h||r, t / [/3(s)]-( d - Q )/ 2 ds 



' J D Q , d 2-( d -«)/ 2 ||7 7 || 1 , t )" n! I n (t,ft) 

= r(i + ft) (A^-w-^iMii.tixi + 

< C* rfJ D(f)"(n!)" ft . 



r(n(l + ft) + 1) 



(The last inequality follows by relation (3.19) of [15].) The first statement 
follows by Lemma [4761 The second statement follows since, 



00 \n 00 \n 00 rxni'-Al 



n=0 n=0 n=0 v ; 

and the last sum is finite for all A > 0, since 1 + h > 0. 
(ii) In this case, 

a n (t,ri) < Cl d J [Q t]2n n"=i v{si,tj)dBdt < Cl d h\\ltt n = C(t) n , 

and 

£[ e ^,^)] < J2 ^a n (t, V )<J2 - AC(t) 

□ 



, = e 

n=0 n=0 



Proposition 4.8 Suppose thatrj : [0,i] 2 — * M+ satisfies the following condition: 
there exist 7 > and 1/2 < H < 1, such that 

n(r,s)< 1 \r~s\ 2H - 2 , Vr,«e[0,t]. (30) 

(i) If f is the Riesz kernel of order a, or f is the Bessel kernel of order 
a < d, and d < 2 + a, then lim e |o L t . e {rf) = L t {rj) exists in L p (fl), for all p > 2, 
and 

supE[exp(XL t . 5 (v))} < C* H<dt J> (\D(t),l- , for all < A < A (t), 
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where C* H d is a constant depending on H, d and a, 

-y ( d-rv\ 2H ( r] _ fy \-[2H-(d-a)/(2)] 



Ao(<) 



and 

' (i-^)"" 1 ^*rWr(i-A)" !ur * 1 - 9B ^/rf = 2 + a 

oo if d < 2 + a 



(ii) If f is the heat kernel of order a or the Poisson kernel of order a, then 
then lim e |o L^ e {rj) = L t (ri) exists in L p (fl), for all p>2, and 

sup E [exp (XL t , e (v))] < exp(AC(t)), for all A > 0, 
where C(t) — C a ^t 2H /a H ■ 

Proof: The proof is similarly to Proposition [3T6J We use the definition ([25]) of 
a n (t,rj), Lemma I3T41 and condition (|30p . 
(i) We have: 

a n (t,r,)<C* H ^ a D(t) n (n!)( d -")/ 2 < oo. 
The first statement follows by Lemma T4. 61 The other statement follows, since 



n=0 ' n=0 ' n=0 * 

If d — a = 2, then the last sum is finite for all < A < A (t) := 1/D(t). If 
d — a < 2, then then last sum is finite for all A > 0. 
(ii) The result follows, since: 

„ n „ n 

a n (t, V ) < C" d / H V(s 3 ,t 3 )dsdt < C n a , dl n / I s ; - t 3 \ 2H ~ 2 dsdt 

)n 
= C(t) n 



and hence 



^ w i<eU»)^ [ ()] 



ft! •« — ' n\ 

n=0 n=0 



□ 

We now introduce the approximation technique of [18] , which will yield si- 
multaneously the existence of the solution of (JTJ) and some representation formu- 
las for the moments of this solution. We review briefly this powerful technique, 
which has been introduced only recently in the literature. The idea is to smooth 
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the noise W, solve the equation driven by the smoothen noise, and then show 
that the solution of the "smoothen" equation converges to the solution of (fTJ). 
For any s,5 > 0, let <ps(t) = 5 — 1 1[ ,5] (*) and 

W e t ;! = [ [ Mt - s) P e(x - y)dW s , v . 

JO JR d 

Note that the noise W £ ' S can be viewed as a "mollification" of W, with rate 5 
in the time variable and rate \fe is the space variable, since 

H>6 = \ V \j) and p e (x) = Jj=p* (j= 

with ip(t) = l [0 ,i](t) and <j>(x) = (2%)- d / 2 e-^ 2 ^ 2 . (Recall that the function , 
defined by u^ s ' (x) — J Rn ip e (x—y)u(y)dy, is a "mollification" of the function u on 
R", if tp e (x) — e~ n ip(x/e) and -0 > is such that J R „ ip(x)dx = 1.) Therefore, 
this approximation procedure can be regarded as a stochastic version of the 
"approximation to the identity" technique, encountered in the PDE literature. 

We consider the following "approximation" of equation (fl}: 

= -Au e ' s +u £ ' S W s ' S , t>0,xeR d (31) 
dt 2 K ' 

u £ ' S x = Uo(x), x e M d . 

We introduce now the rigorous meaning for the solution of (|3"Tj) , which could 
be derived formally from the mild or evolution version of the equation, by ap- 
plying the stochastic Fubini theorem. 

Definition 4.9 An (Tt)t- adapted square-integrable process u — {u £ tx ;(t,x) £ 
M + x R d } is a solution to {57]) if for any (t,x) £ M+ x R d , the process 

Y*f' e > 5 = l [0 ,t](r) \ \ p t -s{x - y)<ps{s - r)p £ (y - z)u e - s y dyds; (r, z) e K+ x 

JO JR d 

exists, is integrable with respect to W , and satisfies 

u e t ' s x =pMx)+ r [ Y r %*> s < 5 5W r , z . 



By 0, the above definition is equivalent to saying that for any (t, x) € K+ x M. d , 
the process Y t,x ' s ' S exists, u t ' x S i 2 (51), is J^-measurable and 

£«'*P) = E(F)p t u (x) + E<y*™*, DF) HV , VF e B 1 ' 2 . (32) 



Before constructing the solution of (|3ip . we mention few words about the 
notation. If X and Y are random variables defined on (fl, T, P), with val- 
ues in arbitrary measurable spaces A", respectively 3^, and h : X x ^ — > K 
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is a measurable function, we define the random variable: E x [h(X, Y)](u>) = 
J x h(x,Y(uj))(P o X~ 1 )(dx). If X and Y are independent, then 

E[E x [h(X,Y)]] = E[h(X,Y)} = E[E[h(X,Y)\X}], (33) 

where E[ ■ ] denotes the expectation with respect to P, and E[ ■ \X] denotes the 
conditional expectation given X. (This result will be used below with X = B 
and Y = W.) 

We have the following result. 

Proposition 4.10 The process u E ' 

S = Wti; (*, x) eR+x R d } defined by: 



u (x + B t )exp( f f A e r < s y B dW r , y -\\\A 



e,<5.Si|2 

\\HV 



(34) 



is a solution of \31\) , where Apy ,B — L ips(t — s — r)p £ (x + B s — y)ds, and 
B = (B t )t>o i s a d- dimensional standard Brownian motion, independent ofW. 

Proof: The argument is similar to the one used in the proof of Proposition 5.2 
of [18]. We include it in response to the referee's suggestion. To simplify the 
notation, we omit writing HV in || • Wnv an( i (•> ')hv- We also omit writing e, 8 
in A E - S > B , i.e. we denote A E - S > B by A B . 

For every ip £ HV, define F v = e w(v)-\W\\ 2 . N ote that 

E(e w ^) = e ll^ll 2 /2 ; Vv eW. (35) 

Since {F v ;tp £ HV} is dense in D 1,2 (see e.g. Lemma 1.1.2 of 28 ), it suffices 
to prove ([SI for F = F v . Define S t , x (<p) = E(u t ] 5 x F v ). Using (O, 

S t , x { V ) = E[E B [u (x + B t )e w{ - AB ^ AB W 2 / 2 ]e w ^-M 2 / 2 ] 

b U,Jt + RA P W{A B + V )-\\A B + v \\*/2 e (A B >V )K 



= E[E B [u (x + B t )e 



= E[E[u Q (x + Bt)e w(A- +v )-\\A- +v f/2 e (A-, v ) 

Let 



h(B, W) = u n (x + B t )e 



Since B and W are independent, E[h(B, W)\B] = f{B), where 

f(b) = E[h(b,W)] =E[u (x + b t )e w( < Ab+ ^ Ab+ ^ 2 / 2 e( Ab ^] 
= u Q (x + ht)e ^^) E[e w^ +v) -\\A^/^ 

= uv{x + b t )e {A "^\ foranyfo=(6 t )t>oeC([0,oo),M d ), 

where C([0, oo),R d ) denotes the space of continuous functions x : [0, oo) — > R , 
and we used (J35j) for the last equality. Hence E[h(B, W)\B] = u (x + B t )e^ B ^ 
and 

StM = E[E[h(B, W)\B]] = E[u (x + B t )e< AB ^]. 
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By the definition of A B and Fubini's theorem, we obtain: 

(A B ,<p) = a H f Al vVr ,. yl \r-T l \ 2H - 2 ]{y~y l )dydy l drdr l 



= [ V £ ' S (t- s,x + B s )ds, 
Jo 

where V £ ' S (t,x) — ((fs(t — -)p e (x — -),<p). Hence: 



S t , x (<p)=E 



u Q (x + B t )e?q> (J V E - S (t- s,x + B s )dsj 



By the Feynman-Kac's formula (see e.g. Theorem 5.7.6 of [19]), {St, x {<p))t,x 
is a solution of the Cauchy problem: 

dSt ^ ] = lAS t>x (<p)+S t , x (<p)V e > s (t,x), t>0,xeR d 
So,x(<p) = u Q (x). 

Hence, 

St.Av) = PyU (x) + / / p t -s(x - y)S s ,y(<p)V € - 5 (s,y)dyds 

JO JR d 

= p y u (x) + a H E [ Y*?' e ' s ip(r' ' ,z')F v \r ~ r'\ 2H - 2 f{z - z')dzdz'drdr' 

= p y u {x)+E{Y t ' x ' £ ' 5 ,DF v ), (36) 

where we used Fubini's theorem for the second equality above and the fact that 
D r ' tZ 'F v = (p(r' , z')F v for the third equality. This concludes the proof of (|3"2"j) 
for F = F v . □ 

Let B % = (Bl) t >o,i > 1 be independent d-dimensional standard Brownian 
motions, independent of W. Suppose that either (fl8|) or (fT9| hold. For any pair 

with i ^ j, let L B ' B be the random variable defined in Lemma l4.61 with 

„ (r8) -f *hV-s\ 2H - 2 if iJ > 1/2 
^ s) -\ \ {r=s} if If = 1/2 

(By Proposition 13.61 a n (t,rf) < oo for all n > 1, and L B B is well-defined.) 

The following result is the main theorem of the present article. 

Theorem 4.11 (i) Suppose that f is the Riesz kernel of order a or the Bessel 
kernel of order a < d, and either H8\) or H9\) holds. Then, for any integer 
k > 2, we have: 

sup E[{u t i) k ] < oo, for allO<t< t (k),x G R d (37) 

e : 5>0 



26 



where 
t (k) = 



k(k - l)D a , d 2- 2H p 2 H T (1 - ^) 



2H 



-1/(2H-1) 



ifd = 2 + a 
ifd<2 + a 



For any < t < to(2) and x G M. d , the limit iit, x '■— hm e jo lim^^o u t 'x exists 
in L 2 {yi), the process u — {ut. x \ {t,x) G [0, to(2)) x R d } is the unique solution 
of {Ip in L 2 (fl), and 

E[u 2 J = E [u (x + Bl)u (x + B 2 ) exp (if" 82 )] := l2 {t, x). 

If x G M d and t < to(M) for some M > 3, then fim £ jo hm^o E\u\' S x ~ u t,x\ P 
/or all 2 < p < M, and for any integer 2 < k < M — 1, 



E[u k t J = E 



l{u (x + Bl)exp\ L t 



B',B j 



3 = 1 



.l<i<j<k 



lk(t,x). (38) 



(ii) Suppose that f is the heat kernel of order a, or the Poisson kernel of 
order a. Then the conclusion same as in part (i) holds, with to(k) = oo for all 
k>2. 

Proof: The argument is similar to the one used in the proof of Theorem 5.3 
of [18]. At the referee's request, we include all the details for the reader's 
convenience. To ease the exposition, we divide the proof in several steps. 
Step 1. We show that for any integer k > 2, 



l[uo(x + Bi)e W \ J2 {A £ ' S ' B \A^ BJ ) HV 



3 = 1 



l<i<j<k 



(39) 



By ([34]) . can be expressed as 



u t i = E B ' 



uq(x + B\) exp 



JR d 



Apf y ^dW r , y --\\A^\\ 2 nv 



for any i = 1, . . . , k. Taking the product over i = 1, . . . , k and using the inde- 
pendence of B 1 , . . . , B k , we obtain that: 



K;f) fc = U EB 



j^B 1 ,...B k 



Uq(x + B\) exp 
Y[u (x + B % t ) exp 



1 



.1=1 



K'!l Bi dW T>v -\\\A^ Bi \\ 2 HV 



JR d 
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Taking the expectation, and using (|33[) with X = {B , . . . , B k ) := B and Y 
W, we get: 



Let 



E 



IJuofx + ^exp (J* J Rd Ap s v Bl dW r , y - \\\A^ Bi f nr ^ \B 



h(B,W)^T[u (x + B t t )exp( [ [ A e r f' Bi dW r<y -h\A 

i=1 \J0 JR d Z 

Then E[h(B,W)\B] = f(B), where 

" k 

Y[u (x + bl)e 



£,5,5*112 

Wnv 



i=l 



f(b) = E[h(b,W)]=E 

= f[u (x + bi)e-^^ b °^/ 2 E 
i=i 

k I k 

= l[u (x + bi)exp[--Y / \\A e ' 5 ' 

s=i y »=i 

= l[u (x + bi)exp \ Y,{A eAh \A e ^)nv 



W(A E - s - b ')-\\A^ s ^\\ 2 nT /2 



;,i)*ii2 



nv i 



for any b = (6\ ... 1 b k ) with b l = [b\) t > Q G C([0, oo), (We used dSSJ) and 
the fact that W^(y +tp) = W(<p) + W(tp) a.s. for any (p,ip £ TCP, which can 
be checked in L 2 (Q), using the fact that W is an isometry between HP and 
L 2 (tt).) Relation (JMJ) follows, since E[{u t S x ) k ] = E[E[h{B ,W)\B]] = E[f(B)]. 
Step 2. We prove that for any (t,x) eR + xR d with t < i (2), 



lim(A e ' ,5 - Bi ) ^ 5 ' B3 ') W p =Lf 2 f\ Ve>0, Vwefi,. 3 , 

<5|.0 ' 

where Oj,j = {cj G f2; B 1 (lS) and B^ (uj) are continuous} (P(Ojj) = 1). 
Let cj G fiij and e > fixed. For any (si,S2) G [0,i] 2 , define 



(40) 



aH So So _ Sl _ r i)^( i - s 2 - r 2 )\n - r 2 \ 2H 2 dr 1 dr 2 if H > 1/2 
J ips(t - si - r)<ps(t - s 2 - r)dr if H = 1/2 

By direct calculation, using Fubini's theorem and the fact that 

[ [ Pe (x+B i Sl -y 1 )p e (x+Bi 2 -y 2 )f(y 1 -y 2 )dy 2 dy 1 = [ P2s (B i Si -Bl 2 -y)f(y)dy, 

JR d JR d JR d 
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(which can be proved by observing that p £l * p £2 — p £l + £2 ), it follows that 



(A 



e,S,B' j^e,S,B 



IHV 



JO 
-B'.B 1 
J t,2e 



Vs(si - s 2 )p2e{B l Sl - Bl 2 - y)f{y)dyds 1 ds 2 
)• (41) 



Note that, for any continuous function g : [0,ty 



lim 

.510 



JO 



m ( Sl ,s 2 Ms uS2 ) dsi d S2 = 7 j° /o l-i - s^ 9 ( Sl ,s 2 )d Sl ds 2 if H > 1/2 

I J 9{s,s)ds if /f = 1/2 



In particular, we consider the (random) function <72e defined by: (?2e(si, S2) = 
J RdP2£ (Bl 1 -Bi 2 ~y)f(y)dy = (p 2e * - Bi), for ( Sl ,s 2 ) e [0,i] 2 . (Note 

that <?2e is continuous by Remark 14.41 ) Then, 



(A 



5,5,3* j^e,S,B : > 



IHV = 



r)s(si, s 2 )g 2£ (si, s 2 )dsids 2 



and 



lim (A 

810 



e,5. B' 4 £,<5,B J \ 

,A )hv 



= L 



Jo 92e{s,s)ds 

B\B j 
t,2e ■ 



ff2e(si,S2)dsldS2 if if > 1/2 
if If = 1/2 



Step 3. We prove that for any (£,#) S M+ x M rf with < < f (2), 

{exp((^ £ ^- B * , A e '^ SJ ) W p)} e , 5 >o is uniformly integrable, (42) 

Suppose first that H = 1/2. Then 775 satisfies condition (|29|) with 1 1 77,5 1 1 1 , t < 1 
(see p. 318 of [H]). By applying Proposition ^. 71 and using (|4ip . it follows that 
if / is the Riesz or Bessel kernel, 



supi? 

£>0 



exp X(A 



ie,i5,B* £e,5,B 



IHV 



<Q d <& AD(i),l- 



VA > 0, 



whereas if / is the heat or Poisson kernel, 



sup-E 

£>0 



exp ^X(A"^ B \A^ B3 )nv)} < e AC(t) , VA > 0. 



Note that both constants D(t) and C(t) depend (linearly) on 1 1 77,5 1 1 1 , t (which is 
bounded by 1), and the function $(2:, a) is increasing in x. We infer that there 
exists an upper bound for the above supremum over e, which does not depend 
on 6. More precisely, denoting by D(t), C(t) the respective constants D(t), C(t), 
in which ||%||i,t is replaced by 1, we infer that if / is the Riesz or the Bessel 
kernel, 



sup E 

£,5>0 



exp 



S,B Z £E,8,B 



IHV 



<C* d $(\D(t),l 



d — a 



VA > 0, 
(43) 
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whereas if / is the heat or the Poisson kernel, 

exp (x(A e ' S ' B \A e ' S ' Bi ) nv )] < e xc{t \ VA > 0. (44) 



sup E 

e,5>0 



Relation ([4"2"]) follows by a well-known criterion (see p. 218 of [5]) by taking an 
arbitrary A > 1. 

Suppose now that H > 1/2. Then r)$ satisfies condition ([50)) (see relation 
(5.13) of [IE])- By applying Proposition 14.81 and using (|41[) . it follows that, if / 
is the Riesz or the Bessel kernel, 



sup 2? 

£>0 



exp 



(A(A 



e,S,B z AeJ^s 



< C*n d<a $ XD(t), 1 - 



d-. 



for all < A < \o(t), whereas if / is the heat or the Poisson kernel, 



supi? 

e>0 



vq>(\{A s > & > B \A e ' s ' Bi )nv 



<e AC «, VA>0. 



The constants D(t) and C(t) depend on 7, which depends only on H. From 
here, we infer that, if / is the Riesz kernel or the Bessel kernel, 



sup E 

e,S>0 



exp (x(A^ B \A^ B3 )nv)} < C* HAa $ (xD(t),l 



d- 



< 00, 
(45) 



for all < A < Ao(t), whereas if / is the heat kernel of the Poisson kernel, 



sup E 

£,<5>0 



exp (\{A e ' 5 ' B \A E ' s > Bi ) nv y\ < e xc ^ < oc, VA > 0. (46) 



Relation (|4*2")) follows as before, noting that 1 < Ao(t) (since t < t (2)). 
Note that ([37)) is obtained as a by-product of ([39)) and ([43)) - ([46)) . since: 



sup E[{u £ t: i) k ] < \\u \t sup E 

e,<5>0 e,S>0 



exp 



HV 



, (47) 



and 



k(k - 1) 



< Xo(t) if and only if t < to(k). 



Moreover, using |[4*3" )) - |[4"T)) , and the fact that D(t),C(t) are increasing func- 
tions of t and Ao(t) is a decreasing function of t, we conclude that, for any 
< T < t (k) and for any (t, x) E [0, T] x M d , 

sup E[(uH) k ] < \\u \LC* H ^ ( k -^±D{T),l- d —^ 

e,S>0 \ Z 1 



if / is the Riesz kernel or the Bessel kernel, and 

sup sup E[(ul' 5 x ) k ] < |K||^exp 

(t,x)G[0,T]xR d s : S>0 



k{k - 1) 



C(T) 
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if / is the heat or the Bessel kernel. Hence, for any < T < to(k), 



sup sup E[(u t f x ) k ] < K T {k) < oo, 

(M)e[0,T]xR d £,5>0 

where Kxik) is a constant depending on uq, k, H, d, a and T. 
Step 4- We prove that for any < t < to(fc), 

limlimM(u^) fc l = j k (t,x). 

First we note that, for any e > 
k 



(48) 



(49) 



lim E[(u £ t i) k ] =E 

8[0 



3=1 



. l<i<3<fc 



:= Ht,x,k{.z)- 



(This follows by Theorem 16.14 of 0, using (35]), gDJ), and @2J|). Next, we 
show that: 

limff Mifc (e) = j k {t,x). 

ej.0 

For this, let (e„) J, be an arbitrary sequence. We have to show that: 

lim Ht,x,k( £ n) = 7k(t,x). (50) 

n — >oo 

We use the fact that in a metric space, x n — > x if and only if for any 
subsequence TV' C N there exists a sub-subsequence N" C N' such that x n — > x, 
as n — > oo, rt e AT". Let JV' C N be an arbitrary subsequence. 

By Lemma [4.61 the limit hf ' B := lim e |o Lf 2 ]f exists in L 2 {VL). (In this 
lemma, we take r/(r,s) = a#|r — s\ 2H ~ 2 if H > 1/2, and T](r,s) — l{ r=s } if 

H = 1/2.) Hence Lf 2 ^ ~ * if ' B m probability, asrn oo,rt 6 AT', and there 
exists a sub-subsequence AT" C A 7 "' such that 

T B\B j T B\B j _ Ar // 

L t — * _L t a.s. as n — > oo, n £ N . 

Note that i?[exp(ALf 2 ^ S )] < oo for all A > 0, respectively for all < A < Ao(i), 
with Ao(<) > 1. (If H = 1/2, we use Proposition 14.71 with rj(r,s) — l{ r=s }. 
If H > 1/2, we use Proposition 14.81 with r)(r,s) — an\r — s\ 2H ~ 2 .) Hence, 
{exp(Lf 2 ' B )} £ >o is uniform integrable. By Theorem 16.14 of [5], 

H t , x ,k( e n) —> lk(t,x), as r»->oo,nSJV". 

Relation (|50[) follows using the above-mentioned subsequence criterion. 
Step 5. We prove that for any < t < to (2), 



lim lim E[ 

E,e' 10 S,S' 10 



e,8 e',5'i 



l2(t,x). 



(51) 
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Similarly to (f3l?P and (f^D|) . one can prove that: 

E[<; 5 Xf] = E[u a (x + Bl)u (x + B 2 t )eM(A eAB \A^^ B2 ) nv )l 

and 

lim (A e ' s ' B \A e '' 5 '' B2 )uv = Lft'£ Ve > 0, Ve' > a.s. 

5,S'\.0 - £+£ 

Relation ([5T|l follows using the same argument as in Step 4 (based on (|4'2p). 
Step 6. We prove that for any t < to (2) and n e f*, the limit 

u t x '■= limhmnf'f exists in L 2 (Q). 

Let t < t (2) and x G K d be fixed. From flSTJ, we obtain that: 

lim lira E\u/ T -u'f\ 2 = 0. (52) 

Let (e n ) n i and (<5„)„ J. be arbitrary. From (|52")) , it follows that (U(™' 5 ")„ 
is a Cauchy sequence in L 2 (VL). Hence, there exists Ut x G L 2 (il) such that 
i?|?Xj™'' 5 " — ut. x \ 2 — * 0. The fact that ut^ does not depend on (e„) n and ((5 n ) is 
proved by a standard argument (see for instance, the proof of Lemma l4~6j) . 

Step 7. We now prove ([38]) . 

If t < t (M) for some M > 3, then sup e s>0 E[(u e t ' x ) M ] < oo, and hence, 
{(M('^) p } £j( 5>o is uniformly integrable, for any 2 < p < M . Since u\' S x — Ut, x — * 
a.s. (along a subsequence), we conclude that -E|?4'f — u t , x \ p — ► 0. In particular, 
lim ei0 hm^o £[K;*) fc ] = F>£j for any 2 < fc < M, and {38]) follows by (09)). 

Step 5. We prove that u = {lit,*; (t, x) G [0, t (2)) x R d } is a solution of (J). 

Let )t,a; be the solution of equation pip (whose existence is guaranteed 
by Proposition I4.10p . Let (e n ) J. and (5 n ) J. be arbitrary. By Step 6, 
— * Ut >x in L 2 (Q), and hence Ut^ 5 ™ — ► u t , x a.s. along a subsequence 
N' C N. Since u\ n x n is JFt-measurable, it follows that u^ x is .^-measurable. 

Let (t, x) G R+ x R d be fixed with t < t (2). We have to show that © holds. 
As in the proof of Proposition ^. 101 it suffices to show that © holds for F = F v 
with ip G TCP. Moreover, it suffices to take ip — lr 0iO ]<^, with a G R+, G V{E. d ), 
since the class of these functions is dense in TiV. 

The idea is to take the limit in (I32p . as e 1 0, <5 J. 0. On the left-hand side, 
E(F v u e t ' x ) — > E(F lp u ttX ), since -E|?4'f — Ut,a;| 2 — > 0. On the right-hand side of 
1(52)1, we have: (see 1(55)1) 

JO JR d 

We have lim ei0 lim 5i0 E(Fu^ s y ) = E{F v u s , v ), for all (s, y) G [0, t] x R d . Note 
that y 6 ' 5 = ip e .s * H, where ip e ,g(8,y) = <ps{s)p E (y) and 

ff (s, y) := a H [ [ <ptf ', y')\s - s'\ 2H *f(y - y')dy'ds'. 

ii + Jm. d 
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Hence, lim^o lim^o V e < s (s, y) = H(s,y), for all (s,y) e [0,t] x R d . Therefore, 
by the dominated convergence theorem (whose application is justified below), 

\\m\im E{Y t ' x ' £ ^,DF v ) HV = f f Pt ^ s (x - y)E(u s , y F v )H(s,y)dyds. 

A direct calculation shows that the limit above coincides with E(Y t - x , DF v )-n-p 
(using the fact that D s ^ y F v = <p(s,y)F v ). 

It remains to justify the application of the dominated convergence theorem. 
Using the Cauchy-Schwartz inequality and (|48|) , we have: 



sup sup \E(u$F v )\ < {K t {2)yf 2 {E{F*)Y/ 2 := K* t . (53) 

(s,Jf)e[0,t]xR d e,S>0 

Note that V £ ' 5 (s,y) = J s (s)I £ (y), where 

Js(s) = ((fs(s- -))l[o,a]>w and I E (y) = (p £ (y ~ ■)> <f>)v(n*y 
Let (3 — (/3t)t>o be a fBm of index H. Then 

Js(s) = =\a H I I \r-r'\ 2H - 2 drdr' = -E[{f3 s - I3 s . s )l3 a ] 

= ±[ S 2H + \ S -6-ar-\ S -ar-\s-6n. 

Using an argument similar to (5.13) of [18], one can show that: 

Js(s) <2Hs 2H - 1 +c a . H , Vs e [0,i],V<5 > 0, (54) 

where c 0) jj = H(a 2H ^ 1 + 1). (This argument is based on treating separately 
the cases a > s, a < s, and considering in each case several intervals for 5.) 
We claim that: 

\I s (y)\<ct, VyGR d ,Ve>0, (55) 

where c$ is a constant depending on d, f and 4>. To see this, we assume without 
loss of generality that <j> = ip * P b, for some ip € C yo (W d ), b > (since functions 
of this form are dense mT(M. d )). Then, assuming that /i(<i£) = g(£,)d£,, we have: 



\ie(y)\ 



(2n)- 2d / e-°M''*W{t)\e- m ' /2 9(W 



< / \T Pe ( y --)(o\ imomn 



vi,) I w(t)\ 2 9(m) 1/2 ( f e-w'gim 1 " 



(2n)-^ 2 U\\ P{Md) [V> (1) (M)] 1/2 
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where we used the Cauchy-Schwartz inequality for the last inequality above. 
Combining (|54|) and (f55|) . we get: 



sup sup \V e >\s, y)\ < c^Hs 2 "- 1 + c aM ). (56) 

(s,y)e[0,t]xR d e,5>0 

From ([53]) and ([5(1, we infer that for any e > 0, S > 0, (s, y) E [0, i] x R d , 

\pt- a (x - y)E(ul : s y F v )V £ - s ( S , y)\ < *(«, y), 

where \Jj(s, y) := K^c^p t - a (x — y)(2Hs 2H ~ 1 + c a _ H )- Finally, we note that * is 
integrable on [0,t] x R d , since: 

/ / t>(s,y)dyds = f {2Hs 2R - 1 + c aM ) ( f Pt-s(x - y)dy) ds 

JO JM d JO \JR d J 

= K* C(t> (t 2H + c aM t) < oo. 

□ 
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A Correction to Theorem 3.13 of [1] 

Theorem 3.13 of pQ gives the necessary and sufficient condition for the exis- 
tence of the solution of the equation: u t = 5 Ait + W in (0, T) x R rf , with 
u(0, •) = 0. This condition is equivalent to saying that ||<7ta;||wp < 00, where 
g tx (s,y) = [2ir(t~ S )}- d / 2 e^p{~\x-y\ 2 /[2{t-s)}} =p t - s (x-y). The condition 
is incorrectly stated in the case of the Bessel kernel, the heat kernel, and the 
Poisson kernel. We state below the correction of this result, whose proof will 
appear as an erratum in [2]. In connection with the present article, we observe 
that: ^ ^ 

\\9tx\\nv=^H f [ \r-s\ 2H - 2 I f (r,s)drd S = ai (t), 
Jo Jo 

since I f (r, s) := / Rd J Rd g tx (s, y)f(y - z)g tx (r, z)dydz = ip*^(r, s). 

Theorem A.l (i) If f is the Riesz kernel of order a, or the Bessel kernel of 
order a, then WgtxWnv < 00 if and only if H > (d — ol)/A. 

(ii) If f is the heat kernel of order a, or the Poisson kernel of order a, then 
WgtxWnv < 00 f or an U H > 1/2 and d > 1. 
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